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Abstract 

We demonstrate how holonomy corrections in loop quantum cosmology (LQC) prevent the Big Rip singu- 
larity by introducing a quadratic modification in terms of the energy density p in the Friedmann equation in the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space-time in a consistent and useful way. In addition, we in- 
vestigate whether other kind of singularities like Type 11,111 and IV singularities survive or are avoided in LQC 
when the universe is filled by a barotropic fluid with the state equation P = —p — /(/?), where P is the pressure 
and f(p) a function of p. It is shown that the Little Rip cosmology does not happen in LQC. Nevertheless, the 
occurrence of the Pseudo-Rip cosmology, in which the phantom universe approaches the de Sitter one asymp- 
totically, is established, and the corresponding example is presented. It is interesting that the disintegration of 
bound structures in the Pseudo-Rip cosmology in LQC always takes more time than that in Einstein cosmology. 
Our investigation on future singularities is generalized to that in modified teleparallel gravity, where LQC and 
Brane Cosmology in the Randall-Sundrum scenario are the best examples. It is remarkable that F(T) gravity 
may lead to all the kinds of future singularities including Little Rip. 

PACS numbers: 98.80.Qc, 04.20.Dw, 04.60.Pp 

Keywords: Loop quantum cosmology; Dark energy; Future singularities 
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I. INTRODUCTION 



The accelerated expansion of the present universe has been supported by various cosmological 
observations, for example, Type la Supernovae DJLD , baryon acoustic oscillations (BAO) yfl, large 
scale structure (LSS) J2I, cosmic microwave background (CMB) radiation J4-6], and effects of weak 
lensing As representative procedures to explain the late time acceleration, the first is to assume 
the existence of an unknown matter called "dark energy" in general relativity, which has negative 
pressure (for recent reviews, e.g., see [|8l-ll2h. The second is to modify gravity, the simplest way of 



which is F(R) gravity (for recent reviews, for example, see Q13HL7I])- 

In general, dark energy can be assumed to be a perfect fluid with the equation of state P = p— /(/?), 
which realizes the current cosmic acceleration. Moreover, the Wilkinson Microwave Anisotropy 
Probe (WMAP) observations indicate that the central value of the equation of state (EoS) is given 
by w = P/p = —1.10 [6D. This means that our universe would be dominated by "phantom energy" 
(/(p) > 0, i.e., w < —1). However, the classical solutions of general relativity for a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model containing dark energy lead, in general, to future singu- 
larities [|18l-l22|l such as the Big Rip and future Sudden singularities. In Ref. [|20j], the finite-time future 
singularities have been classified into four types. 

Recently, the possibility to avoid these future singularities has been studied in the literature, by 



using different approaches to quantum cosmology like loop quantum cosmology (LQC) [|23 



270 (for 



a recent example of a related study on LQC, a super acceleration in LQC and its possible pha se 
transitions, i.e., the crossing of the phantom divide line w — — 1 has been discussed in Ref. 112811) . 



semiclassical gravity, modified gravity, brane cosmology, and etc (for reviews on LQC, see, e.g., [|29 
350). 



Moreover, as an alternative gravitational theory to general relativity, "teleparallelism" Q36 



39] 



has recently been considered. It is constructed by using the Weitzenbock connection. Hence, the 
action is described by the torsion scalar T and not the scalar curvature R as in general relativity 
formulated with the Levi-Civita connection. It is known that the modified teleparallel gravity, so- 
called F(T) gravity, can realize both inflation [40] and the late-time cosmic acceleration |J4 ll— 143|1 1 



1 If there exists a scalar field with a non-minimal coupling to the torsion scalar T even in pure teleparallelism and not in 
F(T) gravity, the late-time cosmic acceleration can be realized [44]. Moreover, the generation of large-scale magnetic 



fields from inflation in the framework of pure teleparallelism has been investigated in Ref. 14511 . 
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(for more detailed references on F(T) gravity, see, e.g., lll2ll ). In Ref. [|46ll . models of F(T) gravity 
in which the finite-time future singularities appear have been reconstructed (a similar investigation 
has been executed also in Ref. |47|]). 

In this paper, we study the features of dark energy cosmologies in the context of LQC. It is shown 
that in conventional cosmology, at the dark energy dominated stage, the domination of dark energy 
finite-time or infinite-time (such as the Little Rip [1481 — iSOM and Pseudo-Rip 115 ill scenarios) future 
singularities may occur, but that in LQC, some of these singularities may be avoided. 

In this theory, holonomy corrections lead to p 2 correction term in the Friedmann equation at high 
energies fl520, which give rise to a bounce when the energy density becomes equal to a critical value 
of the order of the Planck energy density. Thus, this bounce prevents the existence of Rip singular- 
ities in LQC. Indeed, the Friedmann equation in LQC is an ellipse in the plane (H, p). This means 
that the universe moves along this ellipse and the values of the variables H and p are constrained to 
be bounded ll53l - l55n . This is the main difference from the FLRW cosmology, where the Friedmann 
equation describes a parabola in the plane (H, p). Accordingly, in the FLRW cosmology, the universe 
moves along this parabola and the variables H and p can take unbounded values. Hence, for dark en- 
ergy models described with those EoS, Rip singularities could occur. On the other hand, singularities 
that appear for finite values of the energy density and the Hubble parameter, like sudden singularities, 
could still happen in LQC. 

Furthermore, we generalize our analysis on future singularities to F(T) gravity. For this procedure, 
LQC is one of the best examples. We show that in this theory, the universe moves along a curve in 
the plane (H, p), given by the corresponding modified Friedmann equation, and that the dynamics is 
described by the conservation and the corresponding modified Raychauduri equations. Finally, we 
present several non trivial examples of the EoS where future singularities appear. We use units of 



fci 



h — 1 and denote the gravitational constant 87rGN by k 2 = 8n/M P 2 = 1 with the Planck 



mass of Mpi = G 



-1/2 



1.2 x 10 iy GeV. 



The paper is organized as follows. In Sec. II, we review Einstein cosmology (EC) and the singular- 
ities appearing in it. In Sec. Ill, we consider LQC and derive the modified Friedmann and Raychauduri 
equation. In Sec. IV, we study future singularities in LQC and show that the Rip singularities do not 
survive, but Type II and Type IV singularities could still happen. Section V is devoted to investigate 
Pseudo-Rip models, where the universe asymptotically approaches to the de Sitter regime. We illus- 
trate that disintegration of bound structures could occur in LQC, provided that it previously occurs 
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in EC. In Sec. VI, we investigate teleparallel cosmological theories and demonstrate that LQC and 
Brane cosmology in the Randall- Sundrum scenario belong to this category. Finally, conclusions are 
presented in Sec. VII. 



H. EINSTEIN COSMOLOGY 

For the flat FLRW space-time with the metric ds 2 = dt 2 — a 2 (t)dx 2 , where a is the scale factor, 
filled by a perfect fluid with the state equation P = —p — f(p), EC is obtained from the Lagrangian 
£ = \RV — pV where R = 6(H + 2H 2 ) is the scalar curvature, V = a 3 is the volume and H = ^ is 
the Hubble parameter. Here, the dot denotes the time derivative of d/dt. 



This Lagrangian has been constructed in co-moving fluid coordinates (see 115 611 and Sec. Ill C 



of 057Q), and the energy density p has to be understood as a function of the volume V. This relation 
comes from the conservation equation 

d(pV) = -Pd(V) ^p = -3H(p + P)^% = -hp + P) . (2.1) 



dV V 



From P = —p — f(p), we find the differential equation 

dp dV 



(2.2) 



f(p) v 

After integration, we obtain the expression of pas a function of V . For example, if f(p) = —(w + l)p, 
one has 

- y ^ -(1+tU) 



where po is the value of p for V — Vo- 

Note that the Lagrangian can be written as follows: C = — ^yV 2 J — pV. This means that 

the same theory is acquired by removing the total derivative, and eventually we have the Lagrangian 
Ce = — 3V ^ 2 — • T ne conjugate momentum is given by p v = = —^rV = —2H, and thus 
the Hamiltonian becomes 

Re = V Pv -C E = -h V 2 + pV = ~p 2 v V + pV = -3H 2 V + pV . (2.4) 

In general relativity, the Hamiltonian is constrained to be zero. This leads to the Friedmann equa- 
tion 

H 2 = P -. (2.5) 
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This equation represents a parabola in the plane (H, p), that is, the evolution of the universe follows 
this parabola, and its dynamics is given by the so-called Raychauduri and conservation equations. 
The Raychauduri equation is obtained from the Hamilton equation pv = — which reads 

where we have used the conservation and Friedmann equations. Since p v = —2H, we eventually 
obtain the Raychauduri equation H = — |(p + P) = Hence, the dynamics of the universe is 
given by 

2 (2.7) 
P = 3Hf{p), 

provided that the universe moves along the parabola p = 3 if 2 . 

Remark 11.1. Since p = 3H 2 , the dynamics of the system is given by the following 1-dimensional 
first-order differential equations 

H=l^fl. (2.8) 

For the system (12.71) . future singularities are bound to appear in a finite time t s . Actually, they can 
be classified as follows (for details, see [2(1 58]): 



1. Type I (Big Rip): For t — > t s , a — > oo, p — > oo and \P\ —t oo. 

2. Type II (Sudden): For t — > t s , a — > a s , p — > p s and \P\ — > oo. 

3. Type III (Big Freeze): For t — > t s , a — > a s , p — > oo and \P\ — > oo. 

4. Type IV (Generalized Sudden): For t — > t s , a — > a s , p — > and \P\ — > and higher derivatives 
of H diverge. The case that p and/or \P\ tend to finite values in the limit t — > t s is also 



included (590 



There also exist future singularities at infinite time like the Little Rip (LR) Ii48ll . defined as follows: 

H(t) — > oo, when t — > oo. 



The condition that this singularity exists could easily be deduced from Eq. (12.81) and it reads 

f 00 dH 

J H 
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As an example, for the case f(p) = Ap l l 2 , we find []48ll 

H(t) = H e^ A(t - to) , (2.10) 

where -ff = H(t ). 

III. LOOP QUANTUM COSMOLOGY (LQC) 

An approach to quantum cosmology that could avoid singularities is LQC. The main idea is that 
LQC assumes a discrete nature of space which leads, at quantum level, to consider a Hilbert space 
where quantum states are represented by almost periodic functions of the dynamical part of the con- 
nection [30]. Unfortunately, the connection variable does not correspond to a well defined quantum 
operator in this Hilbert space and therefore we re-express the gravitational part of the Hamiltonian in 
terms of almost periodic function. It could be executed from a process of regularization. This new 



regularized 



energies Q52 



iamiltonian introduces a quadratic modification (p ) in the Friedmann equation at high 



6011 . which give rise to a bounce when the energy density becomes equal to a critical 
value of the order of the Planck energy density. 

This theory is constructed as follows. The old quantization of LQC was done by using two ca- 
nonically conjugate variables. One is the dynamical part of the connection, i.e., c, and the other is 
the dynamical part of the triad, namely, p. These variables are related with the scale factor and the 



BOO) 



extrinsic curvature K = |a by the relations (for instance, see 

|p| = a 2 , c = 27A' = 7a, (3.1) 

n 

where 7 = 0.2375 is the Barbero-Immirzi parameter 1)6111 . and their Poisson bracket is given by 
{c,p} = fsgn(p). 

To build the quantum theory in LQC, as a Hilbert space, we usually choose the quotient space as 
the Besicovitch space of almost periodic functions by its subspace of null functions. The Besicovitch 



space (for details, see (16211 ) is the closure of trigonometric polynomials under the semi-norm (in the 
c-representation) 



L 



|*|| 2 = lim / |^(c)| 2 dc 



where c is the connection. Also, all the element of this space have the expansion 
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with fi n G M and a n G ( 2 (the space of square-summable sequences). 

In this space, we can define the operator p as p = — However, the operator c defined by 
c\I/(c) = c\I/(c) is not well defined in this Hilbert space because for a general quantum state ^(c) = 
Y^nez «ne iMnC/2 , we obtain 



|c*l 



L 2 

lim — 



a. 



+00 . 



L— >oo 3 

It is important to stress a key point in LQC. The gravitational part of the Hamiltonian in EC 
contains c. In fact, we have 

3 



Hgrav,EC ~ ~ 3-fiT V — — — C \p 



1/2 



Since the operator c is not well-defined, in order to build the quantum theory, we need to re-define 
the gravitational part of the Hamiltonian. To do that, we introduce the new canonically conjugate 



variables Q34 , 



63J: 



P 



11/2 



IVI 



13/2 



V, 



(3.2) 



with Poisson bracket {(3, V} = |sgn(V). We work in the (3 representation and consider the holo- 



nomies hj(\) = e~ lh ^° ] = cos(^) — iaj sin(^) (see for details Section 2 of ll63ll . Sections II D 
and II E of [|34f| . or Q64Q), where Oj are the Pauli matrices and A is a parameter with the dimension of 
length, which numerical value is determined by invoking the quantum nature of the geometry, that is, 
identifying its square with the minimum eigenvalue of the area operator in LQG. As a consequence, 
it follows A = (see Q63D). 

Thus, using the variables (13.21 ) one has V, gr av,EC = — ^/3 2 V, and since (5 2 , for the same reason 
as c 2 , does not have a well-defined quantum operator, in order to construct a consistent quantum 
Hamiltonian operator, we need an almost periodic function that approaches (3 2 for small values of 



13. This can be executed with the general formulae of LQC fl32i l65l 16611 to acquire the regularized 



671 



Hamiltonian 1164 , 

Hgrav,LQC 



3V 
'^A2 



5>^Tr ^(A)/ ij (A)/ i r 1 (A)/ i - 1 (A)/ ifc (A){^T 1 (A),K}] 



sin 2 (A/3) 



(3.3) 



63 



65 



68]). Note that this 



which captures the underlying loop quantum dynamics (see also 
Hamiltonian admits a quantum version because the operators hj(X) are well-defined in this Hilbert 
space. 
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Then the total effective Hamiltonian is 



H LQC= -ZV^1 + V, (3.4) 



and the Hamiltonian constraint is given by = f. The Hamiltonian equation gives the following 

identity: 

t> fT/o/ i ^dn LQC sin(2A/3) 1 . 

V = {V, 1~Llqcs = ~2 op H = 2lX I 3 = 2~\ arcsm ( 2A 7#)- (3.5) 

By writing this last equation as H 2 = sm a ^ (1 — sin 2 (A/3)) and using the Hamiltonian constraint 
Hlqc = -<=>- s '° 2 ^f^ = f , we find the following modified Friedmann equation in LQC 

H . = £fi_£W_«L + (qit)! = 1 , (3 .6) 

3 V ft/ ft/12 P?/4 
with p c = -^j. This modified Friedmann equation is an ellipse in the plane (if, p) (a parabola in 
EC), and the dynamics is given by the system 



2 V Pc 



P = 3F/(p) , 

where the first equation is the modified Rauchauduri equation in LQC. 



(3.7) 



Remark III.l. The modified Friedmann equation in LQC ( i.e., Eq. A3.6\) ) appears in brane cosmology 



(BC) in the Randall-Sundrum scenario /|52 



3/, 



where the modified Friedmann equation reads 



with A the brane tension. Hence, if we take A = —p c /2, Eq. A3.8\) becomes Eq. A3. (51) . that is, BC with 
negative brane tension is equivalent to LQC. On the other hand, the modified Friedmann equation in 
BC with positive brane tension could be deduced from the Hamiltonian Hbc — ~3V r ^^|^ + V p. 
In the same way, we have obtained Eq. A3. (51) . For the Hamiltonian Hbc> we find 

r2 _ P , P 



H* = r - l + J- 
3 V 2A 



where now A = p c /2. 
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Finally, we remark that for the following parameterization of the ellipse 

H=J^cos( V ), P=^{1 + sin V ), < rj < 2vr , (3.9) 



the dynamics in the case of (13.71) is now given by 

r 1 = x [^F{ri), (3.10) 

V Pc 

where we have introduced the notation F(rj) = /( y (1 + sin 7])). This means that in LQC with a state 
equation of the form P = — p — f(p), we analyze a 1-dimensional first order differential equation. 

IV. FUTURE SINGULARITIES IN LQC 



70HZ6D. 



Future singularities have been studied by using the effective approach of LQC in Refs. (68 
In this section, we examine what kind of future singularities could appear with the effective approach 
of LQC by following our own novel view point. 

A. Rip singularities 

Rip singularities such as the Big Rip and Little Rip are characterized by the fact that H — > oo 
and/or p — > oo. However, in LQC, from the equation of the ellipse in (13.61) . we can deduce that p is 
constrained to belong to the interval [0, p c ] and H in the interval [— \/f§, \/f§]- This means that these 
singularities cannot appear in LQC. For the same reason, the Big Freeze singularity does not occur 
neither. This behavior is very different from those in EC and BC with positive brane tension. Indeed, 
in EC, the Friedmann equation is a parabola p = 3H 2 in the plane (H, p), and in BC with positive 

brane tension (Eq. (13.81 ) with A > 0), the Friedmann equation is the hyperbola jjg = 1. It 

is significant that in both cases, the variables H and p are allowed to diverge, and hence for a model 
with the certain EoS, the Rip singularity appears. In order to clarify this point, we present the helpful 
two examples. 



1. As a first example, we consider the model f(p) = a 2 with a a constant Q69Q. In EC, after 
integration of the system (12.71) . we obtain 

H(t) = y (t - 1 ) + H , p(t) = ~^(t- tof + 3a 2 H (t - t ) + 3H 2 , (4.1) 
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which gives a LR singularity. In BC with positive brane tension (A > 0), the conservation 
equation reads 



which could be integrated and accordingly giving the following result: 



(4.2) 



Pit) = A 




1 + cosh [ 3a y — (t — to) + cosh 



i ( P o - A 
A 



(4.3) 



When t — > oo, we have a LR with the following dynamics: 

A 



H{t) ~ -e 
o 



, p(t) ~ Ae 



(4.4) 



However, if the brane tension is negative (i.e., in LQC), we acquire 



p(t) = A 




1 + sin [ 3a \j —(t - t ) + sin 



i ( P o ~ A 
A 



(4.5) 



that is, p always remains in the interval [0, 2A]. 



2. Another example is given by the model f(p) = Ap x l 2 with A a constant I148h . As seen in 
Eq. (12.101 ), in EC, we have the solution 



Hit) = H e^- t0 \ p(t) = 3H^ Ait ~^ 



(4.6) 



which give rise to a LR singularity. On the other hand, if we investigate the modified Friedmann 
in LQC (namely, Eq. (13.61 )) and substitute H into the conservation equation, in the expanding 



phase (H > 0) we acquire p = VSp^J 1 — y. This equation could be integrated and eventually 
we find 

2" 



P(t) = Pc 



1 



1 _ e Vs\t-t c \ 

I + e V3\t-t c \ 



-oo < t < +oo 



(4.7) 



where t c is the time when the universe arrives at the point (0, p c ). This solution shows that a 
universe moves in an anti-clockwise sense, and bounces (namely, it enters in the contracting 
phase) at time t c . This means that there does not happen the LR in LQC. 
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B. Type II and Type IV singularities 

When f(p) diverges at some point p s , sudden singularities appear because at this point P(p s ) = oo. 
Clearly, this singularity occurs in EC when < p s < oo, and it appears in LQC when < p s < p c . 
In order to compare sudden singularities in EC with those in LQC, we assume that < p s < p c . Note 
that the time when the system arrives at the singularity is different between EC and LQC. Effectively, 
from Eqs. (12.71 ) and (13.71) . we deduce 

f PS d P d P <A 

ts,EC — to+ fo-r/ \ i ts,LQC — to+ . (4.8) 

Jpo VW(p) Jpo yffijl - -e-f(p) 



Equation (14.81) shows that for phantom dark energy (f(p) > 0), we always have t S) EC < ts,LQC 
because . 1 — < 1. 

V P c ~ 

An interesting quantity is the inertial force F inert on a mass m measured by an observer at a 
comoving distance I, and it is given by 

Fi ne rt = ml- = ml(H + H 2 ) . (4.9) 
a 

Near the sudden singularity, since H diverges and H converges, we find 

F inert ^mlH. (4.10) 
Hence, near the sudden singularity, we obtain the following expression in EC and LQC 

Finert,EC = F mert>LQC = TTll^- (l - ^ . (4.11) 

From this last equation, we deduce that for a phantom dark model, F iner ^EC is always positive. How- 
ever, F inertjL Qc is positive for < p < p c /2, whereas it is negative for p c /2 < p < p c . We can 
also conclude that \F inert>L Q \ < F inertjEC , i.e., near the Type II singularity the inertial force in EC is 
always greater than that in LQC. Note that the inertial force is important because when it exceeds the 
(gravitational) force bounding the system then disintegration of bound objects occurs. 

In summary, as shown numerically in early works, Type II and Type IV singularities could survive 
in LQC. Indeed, they survive provided < p s < p c . On the other hand, when p c < p s , the singularity 
disappears and the universe becomes cyclic moving in an anti-clockwise sense along the ellipse. 

Type IV singularities always survive in LQC. This is because they appear at the point (0, 0) and 
near this point the dynamical systems (12.71) and (13.71) coincide with each other. As an example, we 
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consider the case f(p) = — a 2 p 7 with < 7 < 1/4. Near (0,0), the conservation equation reads 
p ~ — \/3a 2 p^ 1 , whose solution is given by 



2 

1-27 



p(f) ^(- ^-^V *.) + ^j ■ (4.12) 

I-27 

This means that the system arrives at the point (0, 0) at t = t s ~ ^(i-^a 2 Po 2 + where all the 
derivatives of higher order than one of H diverge. 



V. PSEUDO-RIP MODELS 



The Pseudo-Rip (PR) has been proposed in Ref. 115 111 . It is defined by P — > —p as t — > 00, 
provided f(p) > 0. In this case, the (mild phantom) universe asymptotically approaches to the de 
Sitter regime. 



As a first example of PR, we consider the model H(t) = H — Hie' 1 * with 7 > 113011 . To obtain 
the EoS corresponding to this example in EC, by inserting this expression into the Friedmann and 
Raychauduri equations and eliminating the variable t, we get 



/(p) = 2flo7(l- Jg^aJ , < p < 3H{ . (5.1) 



P_ 

On the other hand, the corresponding EoS in LQC is obtained in the same way by using the 
corresponding modified equations. The final result becomes 



Pc 



provided 12Hq < p c . 

When PR appears, the quantity of interest is the inertial force F inert . In EC, it is given by 

F inert , EC = ml (M- + P - ) , (5.3) 



and in LQC it is described as 

Since f(p) > (the phantom model), we see that \F inertyL Q C \ < F inertyEC , and thus the disintegration 
of bound structures could occur in LQC, only if previously it occurs in EC. More precisely, let Fq 
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be the force which keeps bounded the structures. Since structures disintegrate when F inert > Fq and 
when p — y pf, we find 

Finert, EC — Tfll—j- , F inertt LQC = ml^- ( 1 — —) . (5.5) 

3 3 V Pc) 

Accordingly, we can conclude that in EC, structures disintegrate if 

F a A. (5.6) 

On the other hand, using the expression is (I3.6I ). in LQC, structures disintegrate for 



provided 

F n <-^. (5.8) 

Note that if the condition (15.81 ) is satisfied, then the condition (15.61 ) will also be met automatically, 
because pf < p c . This proves our statement that disintegration of bound structures could happen in 
LQC, only if previously it occurs in EC. In other words, the disintegration time in LQC is always 
bigger, the corresponding universe is effectively more stable than that in EC. 

Remark V.l. In BC with positive brane tension A, near the PR we get 

Finert,BC = ml^- (l + > F inerttEC , (5.9) 

that is, if disintegration of bound structures occur in EC, then they also happen in BC with positive 
brane tension. 



VI. TELEPARALLEL COSMOLOGICAL THEORIES 

LQC is characterized by two important items: 

1. The canonically conjugate variables (V, (3) are no longer defined by V = a 3 and = . In 
LQC, from Eq. (13.51 ) we conclude that they are defined by V = a 3 and f3 = 4- arcsin(2A7if ), 
and that only when A — y 0, we take the early form. 
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2. From the relation {/3, V} = | -<=>- { V, — = 1, and following the Hamiltonian formulation 
one deduces that V plays the role of the position and — -/3 is its corresponding congugate 
momentum, then Legendre's transformation gives us 

Ulqc = --V/3-£ lqc , (6.1) 

7 

we obtain the Lagrangian in LQC. 

In terms of the variables V and H, it reads 

C-lqc = — — ^- arcsin(2A7iJ) + sin 2 ( - arcsin(2A7# ) J - V> , (6.2) 

which coincides with £ e for small values of A. Note that this Lagrangian does not belong to the 
category of modified gravity models. It is not built from the scalar curvature R = 6 {^H + 2H 2 J 
and the Gauss-Bonnet curvature invariant G = 2AH 2 + H 2 j , i.e., it does not have the form 
£gb = VF(R, G) — V p for some function F. Insteatd of it 



the so-called teleparallel theories of General Relativity [137, 



39|]. 



,QC belongs to the category of 



A. Teleparallelism 

Teleparallel theories are based in the Weitzenbock space-time. To build this space-time, we choose 
a global system of four orthonormal vector fields {e^} related to the vectors {d^} via the relation 
ej = ef In addition, we introduce a co variant derivative V that defines absolute parallelism with 
respect the global basis {e^}, that is, V„ef = 0. From this, we acquire the metric Weitzenbock 
connection T 7 ^ = ejd^e 1 . (Note that this connection is metric and therefore it satisfies V^g^ = 0.) 

The Weitzenbock space-time has identically vanishing curvature (the Riemann tensor vanishes), 
but it is not torsion free. Effectively, we have T 1 ^ u = r 7 „ M — P 7 ^ ^ 0. In order to describe the 
Lagrangian in the Weitzenbock space-time, we introduce the contorsion tensor 

j\ 7 — ^yi j l 7 l j) , 

and the tensor 



Then, we can construct the invariant 
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and define modified teleparallel theories via the Lagrangian of the form C T = VF(T) — Vp. 

The interesting point is tha t if we take the basis {e = d , ei = adi, e 2 = ad 2 , e 3 = ad 3 }, then for 

T = —QH 2 , and hence we acquire Clqc = VF(T) — Vp with the 



the FLRW metric we find Q41 
choice 



F(T) = - 3v arcsin(2A7 V / -r/6) + ^ sin 2 Q a rcsin(2A 7 . (6.3) 

The problematic point in teleparallelism is that it depends on the choice of the global basis. In 
the sense that, if one uses a non-local Lorentz transformation to transform the original global basis in 
another one, in general, one obtains another Weitzenbock connection and thus T could change 11771 
7811 . which does not happen in modified gravity where the invariants do not depend of a global basis. 



However, for a FLRW metric, if we use local Lorentz transformations that only depend on the time, 
that is, of the form = Af (t)ek, then, even though the torsion tensor changes, the torsion scalar 
T remains constant with a value of T = —6H 2 . From our point of view, this gives a consistency 
with teleparallel theories in cosmology (and in particular, with LQC), because T is invariant from 
"isotropic and homogeneous" local Lorentz transformations. 



B. Singularities in teleparallel theories 



Singularities in teleparallel theories have recently been explored in great detail in Ref. [|46Q. Here, 
we present several remarkable results. First, in modified teleparallel theories the modified Friedmann 
equation gives a curve in the plane (H, p) (as we have already seen, in LQC, it is an ellipse, and in 
EC, it is a parabola, whereas in BC with positive brane tension, it is an hyperbola). Thus, if this curve 
is a bounded set, the Rip and Freeze singularities do not survive. 

Now, we consider teleparallel theories with a Lagrangian on the form Ct = VF(T) — Vp. In 
these theories, we have the following modified Friedmann equation 

p = -2F'{T)T + F{T) , (6.4) 

which is a curve in the plane (H,p). Here, the prime denotes the derivative with respect to T. 
Moreover, the dynamics is given by 



H 



f(p) 



< 



4(2F"{T)T+F'(T)) 

(6.5) 



p = 3Hf(p) . 
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This implies that the dynamics of teleparallelism theories are given by a 1 -dimensional first order 
equation. Note that in F(R) gravity theories, the dynamics is more complicated because it is given 
by a first order differential system in IR 3 , i.e., one needs three coordinates, for example, (H, R, p). 
Following this point of view, we investigate a curve of the form p = G(T) for some function G. This 
curve could be obtained from the modified Friedmann equation (16.41) by choosing 



What is interesting here is that for different choices of the curve p = G(T), we find all the kinds 
of singularities. This is contrary to the case of usual LQC which was discussed in previous sections. 
For example, suppose that the function G is monotonic, the dynamical equation (16.51) can be written 
as 

' h = -M(Gr-i)'( p ) 

4 (6.7) 
P = 3Hf(p). 

By taking a dynamic H(t), from the modified Friedmann and Raycahauduri equations 

P = G{T) and H = -^.(GT 1 )' \p) , (6.8) 

we build the corresponding EoS (i.e., /(/>)), which yields the dynamics H(t). As a first example, we 
consider the curve p = T ^ H 2 and the dynamics H(t) = H s e^ H ° { - ts ~ t \ The corresponding EoS is 
given by the function 

In this case, when t — > t s , we find H — > H s , p — > oo and \P\ — > oo, that is, a Big Freeze occurs. 

There another interesting example exists. If we examine the curve p = p s (l — H ^° T ^ and the 
dynamics H{t) = H e Hlt with H , H\ > 0, then the corresponding EoS is given by the function 

/w = ^^0^!, (6 , 0) 

V3H p s 

which leads to a kind of LR, where p does not diverge, because it satisfies p — > p s when t — > oo. This 
is different to the case of third section. 
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VII. CONCLUSIONS 

We have explored the future singularities in LQC. It has been shown that holonomy corrections 
in LQC lead to p 2 correction term in the Friedmann equation and, as shown numerically in early 
works, eventually the Big Rip singularity can be removed. In addition, we have examined whether 
other kinds of future singularities including Type II and Type IV singularities can be cured or not in 
LQC, provided that the universe is filled by a barotropic fluid with the state equation P = —p — f(p). 
Indeed, the Friedmann equation moves in an anti-clockwise sense along an ellipse in the plane (H, p). 
This prevents Rip singularities, but sudden singularizes could survive provided that / diverges at some 
energy density smaller than the critical one. 

Furthermore, by generalizing the above procedures used in LQC, we have studied the future sin- 
gularizes in F(T) gravity. It has been illustrated that from a mathematical viewpoint, this kind of 
theories could be understood as theories modeled by a one dimensional first order differential equa- 
tion. This implies that the study on this theory is easier than that in modified gravity where the 
dynamics of the universe is modeled by a three dimensional first order differential system. 

We also mention that for other models in ordinary F(T) gravity, namely, not in the framework of 
LQC, the Type I and Type IV singularities can eventually appear in the finite time limit for a power-law 
form of F(T). Moreover, the LR and PL scenarios can be realized for specific power-law type models 



of F(T) gravity Il46ll . Accordingly, the features of future singularities occurring in F(T) gravity in 
the context of LQC would be different from those of other models in F(T) gravity. In particular, it is 
very interesting that all the kinds of future singularities as well as Little Rip (infinite-time) cosmology 
are possible in the framework of F(T) gravity. 

Finally, it would be interesting to study the perturbations in LQC by following the recent approach 



developed in Ref. Ir79h . The point is to compare the time required for the disintegration of bound 



objects in the theories with Rip singularity with the time for the future decay of cosmological perturb- 



ations which may be shorter (see feOIP . This will be discussed elsewhere. 
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